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Thermodynamic limit of the canonical partition function with respect to the quark 

number in QCD 

Kenji FukushimE0 

Institute of Physics, University of Tokyo, 3-8-1 Komaba, Meguro-ku, Tokyo 153-8902, Japa^\ 

We investigate QCD in the canonical ensemble with respect to the quark number. We reveal 
that the canonical description in which the quark number is fixed would be reduced to the grand 
canonical description under the thermodynamic limit. Since the grand canonical ensemble contains 
fluctuations of the quark number, the idea of the canonical ensemble is of no use for the purpose of 
defining order parameters for the deconfinement transition. We clarify the origin of such reduction 
and propose an idea to define the order parameter. The Monte-Carlo simulation by means of the 
spin system, which is an effective model of QCD at finite temperature, shows prosperous behavior 
though the results suffer from the severe statistical error due to the sign problem. 

PACS numbers: ll.10.Wx, 12.38.Gc 



I. INTRODUCTION 



Extensive efforts have been shedding light upon the phase structure of Quantum Chromodynamics (QCD) at finite 
temperature, though there remain many subtleties yet. The center symmetry is a prosperous implement to characterize 
the deconfinement transition in pure non-abelian gauge theories (gauge theories without dynamical quarks) The 
order parameter to examine whether the symmetry is broken or not is the expectation value of the Polyakov loop , 
that is, the Wilson line winding around the Euclidean thermal torus. The Polyakov loop vanishes in the confined 
phase, while it takes a finite value in the deconfined phase. In pure gauge theories the critical properties associated 
with the deconfinement transition have been ascertained by the lattice Monte-Carlo simulations. The results are in 
| accord with those anticipated from the center symmetry and the universality argument @]. 
£S| , In contrast to pure gauge theories any definite indicator in order to distinguish the deconfined phase from confined 
one is not established so far yet for systems including dynamical quarks in the fundamental representation. Once the 
thermal excitation of light quarks is allowed, the center symmetry is broken explicitly. If we construct a 3-d effective 
model in terms of the order parameter, namely the Polyakov loop, in the presence of dynamical quarks, those quark 
contributions bring about an external magnetic-like field acting onto the Polyakov loop, which breaks the center 
symmetry p|. Then it is obvious that the Polyakov loop is no longer a proper indicator for confinement because it 
remains finite even in the confined phase as well as in the deconfined phase due to the absence of the center symmetry. 

We note, however, that the notion of confinement should be still articulate even when dynamical quarks are present 
simply because it is an experimental fact. To find out an appropriate indicator to identify the confined phase in 
the presence of dynamical quarks it is essential to clarify where the dynamical screening against the Polyakov loop 
should originate from. In fact, DeTar and McLerran [j| argued that some fractional excitations of dynamical quarks are 
responsible for the explicit breaking of the center symmetry and proposed a new order parameter for the deconfinement 
transition based on the canonical description with respect to the quark triality or quark number. This idea had been 
deserted since Meyer-Ortmanns 6| showed that DeTar-McLerran's order parameter behaves unexpectedly and always 
indicates the confining character for the Z(2) Higgs model on the lattice. Although Meyer-Ortmanns has explicitly 
demonstrated a failure in the proposed order parameter in the model study, it is not apparent in principle what would 
cause DeTar-McLerran's formulation not to fulfill the naive expectation as a proper order parameter. Similar ideas 
for the deconfinement order parameter based on the canonical ensemble have been proposed and discussed by several 
authors also 

0,111 

but those arguments are essentially integrated into DeTar-McLerran's first insight. 
Interestingly enough, we can find a clear-sighted comment in Meyer-Ortmanns' paper as follows: 11 The states at zero 
temperature have a finite particle number hence zero particle density, while the relevant states at finite temperature 
have finite particle density. Therefore one cannot check whether there exist states at zero temperature with fractional 
baryon number by looking at what states contribute to the partition function at finite temperature." In her analysis on 
the Z(2) Higgs model, we can discern no clear embodiment for the statement quoted above. Nevertheless, we would 
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FIG. 1: Typical examples of the thermal excitations of dynamical quarks which respect the center symmetry. The torus 
represents the thermal S 1 and the spatial R 3 . (a) Creation and annihilation process of a quark (q) and an antiquark (q). (b) 
Mesonic excitation composed of qq. (c) Baryonic excitation composed of qqq in the case of iV c = 3. 
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FIG. 2: Typical examples of the thermal excitations of dynamical quarks which break the center symmetry, (a) Single quark 
(q) excitation, (b) qq excitation, (c) Mesonic and single quark (qqq) excitation. 



emphasize that this point is actually the most critical in constructing the order parameter for the deconfinement 
transition. The purpose of the present paper is to give a simple and transparent demonstration to disclose what 
becomes of the partition function with the particle number kept fixed under the thermodynamic limit. 

The problem here can be also stated in more general grounds. The question is whether the grand canonical 
description with zero quark chemical potential is equivalent to the canonical description with zero quark number or 
not. We will find that the canonical description with zero quark number would amount to the canonical description 
with zero quark density whenever the thermodynamic limit is taken. As a result, DeTar-McLerran's idea to project 
out the fractional excitations of dynamical quarks does not work because the states with zero particle density may 
have any fractional excitation of particles which becomes irrelevant eventually in the thermodynamic limit. Then by 
all means can we reach the desired canonical ensemble in which the particle number is fixed at zero? We will propose 
one possibility in the present work. 

This paper is organized as follows: In Sec. [H] we review DeTar-McLerran's idea from the point of view of the 
canonical ensemble with respect to the quark number. Sec. IIIII is devoted to investigating the thermodynamic limit 
of the canonical partition function. Using the simplest example we elucidate analytically and numerically that the 
canonical partition function would be reduced to the conventional grand canonical one under the thermodynamic 
limit. We reveal the cause of the failure in defining an order parameter for the deconfinement transition. In Sec. II VI 
we propose an idea to overcome the problem in connection with taking the thermodynamic limit. Performing the 
Monte-Carlo simulation, we put the idea to the test to confirm that our order parameter would be prosperous. We 
also find that the problem of constructing an order parameter for the deconfinement transition is deeply related to 
the sign problem of the Dirac determinant. The results of the numerical simulations suffer from the severe statistical 
errors due to the sign problem. The concluding remarks are in Sec. ]V\ 



II. IDEA OF THE CANONICAL ENSEMBLE 



A. Canonical ensemble with respect to the quark number 



We will briefly review the idea of the canonical ensemble with emphasis on the possibility to define an order 
parameter for the deconfinement transition. 

First of all, we should notice that the conventional QCD partition function at finite temperature and zero density 
is implicitly formulated in the grand canonical ensemble with zero quark chemical potential, that is equivalent to 
the canonical ensemble with zero quark density under the thermodynamic limit. This fact is easily understood as 
follows. In the imaginary-time formalism of the finite-temperature field theories, the temporal (thermal) extent is 
compactified to the inverse temperature. Then, in addition to the pair excitations of a quark and an antiquark as 
depicted in Fig.n( a ) ; isolated quark excitations as shown in Fig. [3 (a) also satisfy the quark current conservation and 
thus the quark number can fluctuate thermally. In the thermodynamic limit, however, the quark number conservation 
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is recovered because only the states with zero quark number dominate the partition function at vanishing chemical 
potential. This point of view has already been discussed in Ref. an d ^ seems reasonable from the principle of the 
statistical mechanics. 

Usually it does not matter in calculating physical observables whether such fractional excitations of dynamical 
quarks are present. The situation may be totally different for the manifestation of the center symmetry. In the 
language of the lattice gauge theory, the center symmetry is defined as the invariance under the center transformation, 
that is, the transformation from the temporal link variables U<l(x4 = N T ) on a time slice altogether into z-U\{x± = N T ), 
where z is an element of the center of the gauge group (Z(3) for the SU(3) gauge theory). For the quark thermal 
contributions the Z(3) factor z is multiplied every time dynamical quarks wind around the thermal torus. Thus 
thermal excitations gathered typically in Fig.^are center symmetric because z ■ z' = z 3 = 1 for z 6 Z(3). Thermal 
excitations shown in Fig. [21 are, on the other hand, typical examples which break the center symmetry. Then, what 
happens if we exclude any thermal excitation like those in Fig. [2 by using the canonical description where the quark 
number is fixed? Once we employ the canonical ensemble with zero quark number, for example, only such excitations 
as shown in Fig. ^ ( a ) (b) arc allowed so that the canonical partition function would regain the center symmetry. 
Certainly the symmetry itself is explicitly broken in the Lagrangian and nevertheless the physical states remain center 
symmetric. Owing to the recovered symmetry, we can presume that the Polyakov loop will provide a criterion for the 
deconfinement transition in the same way as in the case of pure gauge theories. Essentially, this is the idea proposed 
originally by DeTar and McLerran in Ref. pg. 

The canonical description of QCD at finite density is available by the Legendre transformation from the quark 
chemical potential \i to the quark number n q (In this paper we simply write the quark number to mean the 

quark number minus the antiquark number.) For later convenience, we explain the transformation procedures in some 
mathematical details. To make our calculation tangible, we adopt the Wilson fermion on the lattice from now on. 
The final results, of course, can be understood not only in a specific formalism but in a general way also. 

The quark number density is the fourth component of the conserved current given by 

H(x) = k{$(x){\ - jipiix^ix + 4) - $(x + 4)(1 + 14 )UI (x)i/j(x)}, (1) 

where k is the hopping parameter and the Wilson parameter is chosen as r = 1 for the sake of simplicity. By imposing 
the constraint onto the configurations we can write the canonical partition function with respect to the quark number 
n q as 

Z CE (n q ) = J VUV^e- s ^- s ^ u ^U^^2j 4 (x)~n^j 

^ e -i»«* / 'upp^e-^-*^*'^^, (2) 

where 5g[U"] is the gluon action and Sf[U, V>, ip] is the action of the Wilson fermion with r = 1. The lattice spacing a 
and the number of lattices along the temporal direction N T give the inverse temperature, i.e., (3 = N T a. The quark 
parts of the action are put together as 

X 

X j 



(l+^)^{x)(l-^)U 4 (x)^{x+4) + (l-i^)^(x+4)(l+74)t/i(x)V>(x) 



(3) 



Here i(f> can be regarded as an imaginary chemical potential and the integration over <f> corresponds to the Legendre 
transformation from the chemical potential to the quark number j runs from 1 to 3 in the s pat ial directions. In 
the same way as the usual prescription to treat the chemical potential in the lattice gauge theories [12| we exponentiate 
the imaginary chemical potential as follows; 

0a (.<j>a\ .(pa ( ,cj>a\ 

l + ^exp^-j, l^-^exp (4) 

which makes only higher order corrections of 0(a 2 ). The fundamental reason to demand the above prescription is 
because of the gauge invariance, or the center symmetry in the present case. With these alterations the partition 
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function given by Eq. J5J becomes exactly center symmetric for n q — 0, ±3, ±6, ... (0 mod 3) as it should be. In 
order to see it apparently, transforming the quark fields by 

ift'( x ) = (* i r x *il>(x), $(x) = e-^^ix), (5) 
we rewrite the quark action in terms of the transformed fields as 

-Sf[E^',V'] 

+ ^ [e^^Xl^)^,^^ (6) 

X 

§f[U, '0, if>] is the rest of the action from which the terms involving U^x, N T ) are subtracted. It is obvious that the 
center transformation U&(x,N T ) — » e l2 ^ k ^Ui(x, N T ) is equivalent to the shift of (ft by 2nk/3. We can immediately 
recognize from Eq. (J2J that the canonical partition function is certainly invariant under (ft — > (ft + 2nk/3 as long as the 
quark number satisfies n q — mod 3. 

The integration of the partition function in terms of (ft exactly corresponds to the decomposition of the Dirac 
determinant into each part with fixed quark winding number. This is seen transparently if the Dirac determinant 
obtained from the quark action Eq. ©, that is denoted by det D[U, (ft) from now on, is calculated by means of the 
hopping parameter expansion. The expanded terms are represented by closed quark paths similarly as depicted in 
Figs, ^and |21 Every time the quark path winds around the thermal torus, it picks up the factor e 1 ^ so that the Dirac 
determinant can be decomposed according to the quark number n q \ 

00 

detD[U,(p}= V e in «* det [17], (7) 



71 a — — 00 



where det D^ n ^\U] is the part of the Dirac determinant with the quark path winding n q times around the thermal 
torus. Thus the integration over (ft in Eq. singles out some specific combinations of the operators coming from the 
Dirac determinant under the condition that the quark winding number is fixed. This result is not particular to the 
case of the Wilson fermion formalism on the lattice. In general, the canonical partition function is given by 

Z CE (n q ) = j* T^e-'™^ J T>U det D[U,<f>]e- Sa M ( 8 ) 

with the Dirac determinant det D[U, (ft] under the twisted boundary conditions, 

if)(x,x A = P) = -e^VC^ ^4 = 0), Tp(x,X4, =13) = -e -i *^(»,x 4 = 0). (9) 

It leads to the canonical partition function in the same sense as in Ref. ;8], that is, the partition function with the 
Dirac determinant replaced by its specific part with some fixed quark excitations. It is interesting that the constraint 
onto the configurations, imposed by Dirac's delta function, results in the restriction onto the combination of the 
operators after all. This is because of the non-local character of the quark (fermion) fields. 



B. DeTar— McLerran's order parameter 

Here we will explicitly construct DeTar-McLerran's order parameter for the deconfinement transition. They intro- 
duced the quark triality t q defined as n q mod 3 and considered the canonical partition function with respect to t q . 
It is available from the canonical partition function with respect to n q as 

00 

ZcE,(t q ) = ^2 Z CE (n q = 3m + t q ) 

m— — OO 

= - J2 e ~ iM [l>UdetD[U,<p]e- SG M. (10) 

0=0, ±27r/3 

Then DeTar-McLerran's order parameter is given by 



v M = i^h^i (id 

Et„ Z CE(t q ) 
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for arbitrary Af 7^ 0. When the quark number is fixed instead of the quark triality, we can define also in the same 
way as 



ZcE(n q = Af) 

T,n a Z CE(n q )' 



(12) 



which is expected to behave similarly to Tj^. In the confined phase where all the physical excitations are color-neutral, 
any thermal excitation of isolated quarks would cost infinitely large energy. Therefore it follows that Zce (tq = Af) = 
ZcE(n q = Af) = for Af ^ mod 3. After the dynamical breaking of the center symmetry occurs, fractional quark 
excitations are allowed at finite energies so that the canonical partition function can take a non-vanishing value. 
Hence, Vj\f and Vj\f are anticipated to serve as order parameters for the deconfinement transition. 

We can manipulate the order parameter in other forms comparable with the conventional definition, namely the 
Polyakov loop. As discussed in Refs. [1,0, the expectation value of the Polyakov loop itself can be expected to be a 
proper order parameter in the canonical ensemble. Because the center symmetry becomes manifest in the canonical 
partition function with the quark number n q = mod 3 (or the quark triality t q = 0), the deconfinement transition 
could be characterized by the dynamical breaking of the center symmetry. Thus another candidate for an order 
parameter for the deconfinement transition is given by (trL(x))cE with t q — 0, that is, 



(trL(x)) C E = 



\Y,<t>I vu det D[U, (t>] e- SG M tiL(x) 



Z C E{t q = 0) 



(13) 



In the same way as in the case of Vj\f a similar variation is also possible for (trL(x))cEj *-e., we can consider (trX(x))cE 
with n q = for example as 



(tri(x)) C E = 



# J T>U det D[U, </>] e- s °M tri(x) 



ZcE(n q = 0) 



where L(x) stands for the Polyakov loop defined by 

L{x) — Texp 
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dx4 Ai(xi) 



(14) 



(15) 



Here T represents the time-ordering and A^x) is the temporal component of gauge fields. 

We can explicitly evaluate the integration with respect to <f> once the actual form of det D[U, <j>\ is given. Now we 
will exploit the hopping parameter expansion. Following the calculation of Ref. [l3l | the Dirac determinant amended 
by (j) can be readily expanded in the lowest order of the hopping parameter expansion as follows; 



det D[U,4>] = cxp 



fl"53(e*trL(5) +e- i< ^trL t (x)} 



(16) 



where H = 2(2k) Nt and L(x) here is the Polyakov loop defined on the lattice by 

L(x) = Y[ Ua{x,xa), 

X4—a 



(17) 



which can be graphically represented like in Fig. [2 (a). With the notations M = trL(x) — \M \ e l6M , the integration 
in terms of <j> in Eq. results in [l^ 



Jo 27T n " 



Ml 



(18) 



where I v (x) is the modified Bessel function of the first kind. By differentiating the integrand with respect to </>, we 
can have the expression for the expectation value of the quark number as 



(19) 



This relation can be confirmed analytically also by the formula /„_i(x) — I„ + i(x) = (2v / x)I v (x). Also it can be 
expressed as n q = 2iHlm(J2 tiL(x)e 1 ^) rlq . This form of the relation has been found already in Ref. [l4j. It would be 
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FIG. 3: The schematic picture for the emergence of fractional quark excitations allowed even in the canonical ensemble under 
the thermodynamic limit. 



worth mentioning that a different but similar relation has been discussed also in a numerical way in Ref. |15| . which 
suggests that the above relation could persist beyound the leading order of the hopping parameter expansion. 

In the lowest order of the hopping parameter expansion, we can write down DeTar-McLerran's order parameter as 

^ JVUe~ s -M+^lMmM\) 

14 jT)lf e -S G [U]+2H\M\cos6 M ■ ^ U > 

As for the expectation value of the Polyakov loop in the canonical ensemble, we can immediately write down as 

1 JVUe- s ^+^'I (2H\M\)\M\ 
(tri)cE -y JVUe-Solu] Io{2 H\M\) ' (21) 

In the next section, we address the problem of taking the thermodynamic limit for these would-be order parameters 
(I2U|) and 1)21(1. The point will turn out that the sum of the Polyakov loop M = trL(x) is an extensive quantity and 
as large as proportional to the volume. 



III. THERMODYNAMIC LIMIT 
A. Failure as order parameters 



We consider the thermodynamic limit here, that is, we make the system volume go to infinity. Then we will find 
that both V_\f and (trL)cE fail serving as order parameters for the deconfinement transition. Roughly speaking, the 
failure as order parameters stems from fractional quark excitations permitted even in the canonical ensemble under 
the infinite volume limit. As schematically depicted in Fig. [31 a pair of a quark and an antiquark can be excited 
(drawn by the thin curves) because such an excitation is mesonic and the net quark number is zero. When the distance 
between source quarks (drawn by the thick curves) gets large, it would be regarded as the system composed of two 
mesonic excitations as shown in the left figure of Fig. |3J In the limit of the infinite distance between source quarks, 
two mesonic excitations should become individual (clustering decomposing property as shown in the right figure of 
Fig. |3J) , which means that the Polyakov loop associated with a single quark excitation can be effectively screened by 
the fractional excitation of dynamical quarks. Although any fractional excitation of dynamical quarks is prohibited 
at first by the definition of the canonical ensemble, it would be allowed eventually in the infinitely large volume limit 
due to the clustering decomposing property. As a result, Pjvyo amounts to non-zero V_\f=o and (trL)cE is reduced to 
the Polyakov loop expectation value in the conventional grand canonical ensemble under the thermodynamic limit. 

Let us confirm the above intuitive understanding in the actual expressions given by Eqs. (|2U[1 and l|21|). Without any 
spontaneous breaking of the center symmetry, "Pjv is expected to vanish for M ^= mod 3. This expected property 
derives from the integration with respect to 8m in Eqs. I|20() and l|21l) . Owing to the projection by the integration 
over 4>, the partition function depends on 8m only through the factor e lA ^ M . Therefore, if N takes a non-zero value, 
the integration with respect to 6m makes zero as a result of the whole average over the phase factor. Unless external 
perturbation is introduced, Pa/Vo = should be maintained even in the thermodynamic limit. In other words, any 
spontaneous breaking of the symmetry cannot be described without external perturbation. Thus even after taking 
the thermodynamic limit, we have vanishing "Pvyo m ^ ne absence of external perturbation regardless of whether the 
system lies in the confined or deconfincd phase. This is what Meyer-Ortmanns found in her analysis. 

In the presence of small perturbation denoted by h, the phase factor is modified to c lJ ^ e M+2h\M\cos8 M _ j n ^he 
thermodynamic limit, \M\ may become as large as proportional to the volume and then only the stationary point 
9m = dominates the integration with respect to the Polyakov loop. Consequently loses any M dependence to 
lead to Vj\f ^ for arbitrary Af. In other words, the regained center symmetry is always broken and the partition 
function is reduced into that in the grand canonical ensemble in which the center symmetry is explicitly broken 
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whenever small perturbation is introduced in order to realize the spontaneous symmetry breaking. Thus the situation 
is sharply contrast to Meyer-Ortmanns' analysis. Once external perturbation is applied, the canonical description 
resolves itself into the grand canonical description. If we deal carefully with the thermodynamic limit and try to 
retain the canonical description, even the spontaneous breaking of the center symmetry becomes unavailable. In any 
case, the state of affairs is somewhat of antinomy. 

In order to gain a deeper insight, we will deal with a toy model in which the Polyakov loop dynamics is emulated 
by the spin system in the next subsection. 



B. A toy model 

In this subsection we investigate the thermodynamic limit of the canonical ensemble by means of the Ising model 
as a toy tool. The argument of the center symmetry tells us that the Polyakov loop dynamics can be effectively 
described by the Ising model with the Z(2) symmetry for the SU(2) gauge theory and by the Potts model with the 
Z(3) symmetry for the SU(3) gauge theory. In fact, such an effective description has been established also in the 
numerical study [l6| . 

For simplicity, we will focus here on the Ising model with the Z(2) symmetry, though the discussion is easily 
extended to the case of the Potts model. The counterpart for the order parameter H20|) can be immediately expressed 
in terms of the spin variables as 

E, s) e-'E '»'»i(e* £ •» - e~ H ^ •») 
1 -i ~ y. . e - J £ »»«•»+*£•» { ' 

with the exchange interaction J between the nearest neighbor sites under a magnetic field H. Also the counterpart 
for the order parameter (|21|l can be written as 

; ^_l E,„e-^-S(e^- + e-^..) E ^ 

Woe ~y £ {a) e-JE.»'»!(e*E^+e-ffi:»») ' [ ' 

For the purpose of investigating what happens in the thermodynamic limit, it will be sufficient to consider the 
simplest case in which the exchange interaction is absent (J = 0). Then the order parameters given by Eqs. (|22|) and 
(12311 can be exactly evaluated, namely, they are trivially zero. The naive expectation is that these order parameters 
remain zero while the magnitude of the exchange interaction is small, regardless of the presence of an external magnetic 
field H ^ 0. They will come to take a finite value when the strength of the interaction grows large enough to bring 
about the dynamical symmetry breaking. Actually, however, this naive expectation is not realized at all. We must 
follow the standard procedure to describe the spontaneous symmetry breaking; we add an infinitesimal external field 
denoted by h and then take the thermodynamic limit (V — > oo) and finally turn off the external field (h — > 0). The 
order parameters under an infinitesimal external field h is calculated for J = as 



2 

We can draw the similar result for another candidate, that is 



£ w i( e (W)E'.- e -(")5;'») 

(2cosh(iJ + h)) V - (2cosh(ff- h)) V 
2(2cosh(# + h)) V 

' (V^oo, H>/i>0). (24) 



(s>ct g 



l Ew i( e (W)E H e-'")^-)E^ 

(2 cosh(# + h)) V ~ 1 2 smh(H + h) - (2 cosh(F - h)) V ~ 1 2 sinh(iJ - h) 
(2 cosh(# + h)) V + (2 cosh(iJ - h)) V 
->tanh(i? + h) (V->oo, H^h>0). (25) 

The above results of Eqs. I|24|l and (|25|l clearly signify that the would-be order parameters given by Eqs. (|22|l and 
(12311 are to be disqualified under the thermodynamic limit because they always remain finite. In other words, any 
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expectation value calculated in the canonical ensemble in which the center symmetry is seemingly recovered would be 
reduced into that calculated in the grand canonical ensemble in which the symmetry is explicitly broken. The point 
is that the contribution among the projective superposition with respect to <f> survives only when it is parallel to the 
direction of the external magnetic field. As a result, the projection into the canonical ensemble is diminished only to 
result in the same description as in the grand canonical ensemble. We can intuitively understand the situation in the 
following way: In the thermodynamic limit, the projective superposition in Eq. l)2.Sll for example becomes 

e ff £»»+ e - ff £*» ->e H l£' , *l ) (26) 

which follows that the spin configurations for which e H ^ s " dominates should be definitely separated from the spin 
configuration for which e~ H ^ Sn does, that is, the system becomes non-ergodic. Which term survives in the thermo- 
dynamic limit depends on which direction an external field is introduced in. 

The same mechanism makes the canonical ensemble in QCD be unstable under the thermodynamic limit. For the 
configuration of the Polyakov loop with a macroscopic order, M = ^ trL(x) = \M\ e 10M ~ 0(V) — > oo, the integration 
with respect to <j> in Eq. (|18|l becomes 

p -iig0 . „H(e 1 *M+e - '*M t ) _ p -iri 9 . p 2H\M\ cos(0+0 A /) 



e • e v 1 = / 7, — e 

o 2-7T J 2tt 



27J- Q-m q (j>o . q2H\M\ cos(<Po+8m) 



2w\-2H\M \ cos(<j)o + 6m)\ 



(|M| ~ 0(V) -> oo). (27) 

0o=— Su 



In the last line of the above equations only the leading order of the saddle-point approximation (2H cos(<fo + &m) is 
maximized at 0o = —&m) is left, which becomes exact when \M\ goes to infinity. The meaning of the saddle-point 
approximation in Eq. (|27|l is absolutely the same as what Eq. means. Then infinitesimal perturbation like the 
h field in Eqs. (|24|l and l|25|l forces only specific 6m to be favored in the thermodynamic limit. Thus neither DeTar- 
McLerran's order parameter nor the Polyakov loop expectation value in the canonical ensemble would serve as a 
proper order parameter for the deconfincmcnt transition. 



IV. IDEA TO OVERCOME THE PROBLEM 



A. Prospect as order parameters 

So far we have clarified why DeTar-McLerran's order parameter does not work as expected under the thermo- 
dynamic limit. Nevertheless, we have prospects of making use of DeTar-McLerran's order parameter to identify 
the deconfinement transition. The idea is quite simple. The spin-system counterpart for DeTar-McLerran's order 
parameter given in Eq. (1221 can be equivalently expressed in the following form: 

-plsing _ 2 \° lj_ / q\ 

- — w^Tj — ' ( } 

where (•■■),/ denotes the ensemble average with the action —J^2 s n s m . Let us imagine that we compute the order 
parameter of Eq. (|28|l by the method of the Monte-Carlo simulation. First we generate the ensemble of configurations 
with the probability specified by the action — J^2 s n s m-> which has both the symmetric phase with random alignments 
and the broken phase with a spontaneous magnetization. Next we calculate the ensemble average with those generated 
configurations. In the symmetric phase, the numerator of Eq. I|28|l should vanish, while it take a finite value in the 
broken phase. Hence it is apparent that the would-be order parameter given in Eq. (|28|l will serve as a proper order 
parameter. If we prefer the spin expectation value to DeTar-McLerran's formulation, it would be enough to rewrite 
Eq. (H into 

WOE " V ( e HEs n+e ~HEs n)j ■ ^> 

What is important here is that the probability for configurations is specified by the symmetric part of the action. If 
we adopt the whole action, for example — JJ2 s n s m + lncosh(i? s «) f° r Eq. (J2SJ, to generate the configurations, 
the ensemble of configurations becomes non-ergodic in the thermodynamic limit and the would-be order parameters 
should inevitably fail as demonstrated in Eqs. (|24|l and l|25|) . In other words, the ensemble of configurations given by 
— J^2 s nS m + hi cosh(_ff ^2 s n ) is unstable against configurations with ^2 s n ~ 0(V) because of the long-ranged nature 
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of the interaction term lncosh(iJ s n)- As a result the spin configurations with s n > is completely decoupled 
from those with ^ s n < and the canonical description is reduced to the grand canonical one. 

The above idea is quite simple in itself. In the rest of this subsection, we will discuss the physical meaning of the 
above prescription to generate the ensemble of configurations with the probability specified by the symmetric part of 
the action. 

Looking back at the argument around Eq. (|27(l . let us reconsider the validity for the saddle-point approximation 
which is expected to be exact in the thermodynamic limit. Once \M\ is really as large as of order 0(V), the argument 
given in Eq. (|27[1 has no suspicious controversy at all. However, going back to the second line of Eq. we realize 
that the coefficient in front of <fi is at most n q which is now of order 0(1). If we considered the thermodynamic limit 
where we take n q ~ 0(V) — > oo with the number density n q /V fixed, the aforementioned argument of the saddle-point 
approximation would become exact and the calculation of the canonical ensemble is absolutely identical as that in 
the grand canonical one, as closely discussed in Ref. ^t|- Then the expectation value of the phase of the Polyakov 
loop is finite and almost proportional to the quark number density [l5}, as seen in Eq. (|19fl . 

As we have already emphasized, the canonical ensemble for the purpose of characterizing the deconfinement tran- 
sition is different. What we should be confronted with is the ensemble where n q is fixed at so small number of order 
0(1) that the saddle-point approximation would be inapplicable. So as to realize this situation an additional condition 
is needed; \M\ — \ ~^2^L(x)\ must be kept of order 0(1). This condition can be achieved in the form of the ensemble 
average in Eqs. (|28|) and l)29|) as numerically demonstrated in the next subsection. 

We shall summarize our findings and assertions here again. 

• Even when the quark number is fixed in the canonical formulation, it is the quark number density that is 
actually fixed under the thermodynamic limit because the system becomes unstable against configurations with 
a macroscopic order. Thus any fluctuation of the quark number of order 0(1) is allowed even in the canonical 
description, though it is unintended. As a result the notion of confinement becomes obscured. 

• The idea is simple: If we completely get rid of the fluctuation of the quark number, which is negligible in the 
thermodynamic limit but responsible for screening the Polyakov loop, the canonical ensemble with respect to 
the quark number is expected to recover its meaning. Then DeTar-McLerran's order parameter and also the 
Polyakov loop should serve as order parameters for the deconfinement transition. 

• The fluctuation of the quark number is induced by the infinitely non-local interaction arising from the projective 
superposition. Such fluctuation can be excluded by taking the ensemble average with configurations generated 
by the symmetric part of the action. The asymmetric part of the action is regarded as included in the operator 
part to be averaged over. 

B. Numerical tests 

By using the Monte-Carlo simulation, we numerically calculate Eq. (|29[1 which can be comparable with an ordinary 
magnetization in the presence of an external magnetic field. The lattice size is chosen as 100 x 100 and we adopt the 
standard Metropolis algorithm. The strengths of the exchange interaction and the magnetic field are J = 0.5/T and 
H = 0.1/T respectively, where T is the temperature. The thermalization of the spin system is well achieved after 
10000 times sweeps, as can be confirmed by the magnetization distributions shown in Fig. 0] for T = 0.5 and T = 2.0. 
It should be noted that the magnetization is typically of order 0(V) in the ordered phase (T = 0.5) and of order 0(1) 
in the disordered phase (T — 2.0). 

The resultant (s)^ 8 is presented in Fig. [S] Filled circles stand for the spontaneous magnetization in the case 
of H = and indicate that the spin system goes through the second-order phase transition at T c = 2.27J = 1.13 
(theoretical value). Blank circles represent the magnetization in the presence of the magnetic field H ^ and 
show a smooth crossover as the temperature raises. We plot (s)c]^ s by crosses with error-bars. Although we take 
10000 samples by each 10 intervals, the statistical error is still large because the term e H ^ s ™ + e~~ H £ s " in both the 
numerator and the denominator of Eq. (|29|) can become so large that the operator to be taken the ensemble average of 
must furiously fluctuate to lead to large statistical errors. In a sense we can say that this problem of slow convergence 
arises from the same origin as that of the sign problem for simulations at finite density, in particular, the situation is 
alike the Glasgow method in lattice QCD at finite density. This parallelism to the sign problem will be more apparent 
for the order parameter given by Eq. I|28|) or Eq. (|20l) . As discussed in Sec. IHII the canonical partition function can 
vanish due to the average over the phase factor coming from the Dirac determinant. 

The results shown in Fig. [5] seem to suggest that the would-be order parameter (s)ce S should work as expected. 
Unfortunately, at this stage, we cannot address any stronger statement due to large statistical errors. In order to 
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FIG. 4: Magnetization at each sweep step after the thermalization achieved by 10000 times sweeps. The left figure is for T = 0.5 
in the ordered phase and the right for T = 2.0 in the disordered phase. 
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FIG. 5: Behaviors of the magnetizations and the order parameter as functions of the temperature. 

improve the accuracy, some essential ingenuity such as a reweighting method 0] are needed rather than simple 
accumulation of many samples. That is beyond the scope of the present paper. 



C. Construction of the order parameter in QCD 

It is easy to write down the QCD counterparts for the order parameters given by Eqs. I|28|l and (|29|) in the spin 
system. Here we present the expressions only for DeTar-McLerran's order parameter defined by Eq. (|llfl and the 
Polyakov loop expectation value given by Eq. I|13fl . The formulation can be straightforwardly extended to general 
cases. 

The quark contribution from the Dirac determinant can be expressed as the effective action in the following form; 

det D[U] = e - s ^ u ' u ^ Xi=N ^\ (30) 
where the link variable Ui(x4 = N T ) receiving a modification under the center transformation is explicitly written for 
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convenience. We can decompose the effective action according to the triality, that is, 

S^[U] = l X) e- a * kn / 3 S eS [U,zZU 4 (x 4 = N r )] (k = 0,±l) (31) 
n=0,±l 

with zo = e l27r / 3 G Z(3). Among the decomposed parts, all the contributions with non-trivial triality (k = ±1) are 
bound to be included in the operator to be averaged over. Therefore the would-be order parameter given by Eq. 111|) 
can be written for Af — 1 as 

1 / CjC 1 ) of" 1 ) _1 . c(!) .-lot-l) ,-lc(l) , 0<-!)\ 

i/e — o* — o ff + Zq e~ =ff ~ o «« + zq e~ ° <=« 0,5 cff \ 

A = 7 ^nrs — ■ ( 32 ) 



sym 



The ensemble average is to be taken with the probability weight specified by the symmetric part of the action, i.e., 
—Sg[U] — S^[U]. In the same way, the expression for Eq. I|13[) can be written as 

( (e rft °cft -(- c z °°cff z o °cff -(- c o °ctf z o°otf ) trLfa;) ) 
(trL(f)) CE = 22. . (33) 

/ sym 



V. CONCLUDING REMARKS 



In order to identify the deconfinement phase transition we employed and elaborated the idea of the canonical 
ensemble with respect to the quark number or triality. We clarify why the canonical ensemble would be reduced 
into the grand canonical ensemble eventually in the thermodynamic limit. In order to overcome the problem of the 
thermodynamic limit, we propose a prescription to compute the ensemble average. The idea is tested by means of 
an effective model in terms of spin variables in the presence of the magnetic field. The results seem prosperous but 
turn out to suffer from the severe sign problem. Although our definition of the order parameter for the deconfinement 
transition can be applied to the Monte-Carlo simulation on the lattice, it is inevitably necessary to get over the sign 
problem inherent in the lattice QCD at finite density. We find it interesting that two serious questions - one is the 
criterion for the deconfinement transition at finite temperature, and the other is the sign problem in the lattice QCD 
at finite density - are closely related to each others. 
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